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A general reformulation of the exchange energy of 5/-shell is applied in the analysis of the magnetic 
structure of various actinides compounds in the framework of LDA+U method. The calculations 
are performed in a convenient scheme with essentially only one free parameter, the screening length. 
. . . The results are analyzed in terms of different polarization channels, due to different multipoles. 

0^ I Generally it is found that the spin-orbital polarization is dominating. This can be viewed as a 

. strong enhancement of the spin-orbit coupling in these systems. This leads to a drastic decrease in 

' spin polarization, in accordance with experiments. The calculations are able to correctly differentiate 

magnetic and non-magnetic Pu system. Finally, in all magnetic systems a new multipolar order is 
f— I , observed, whose polarization energy is often larger in magnitude than one of spin polarization. 
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The magnetism of actinide systems shows a very rich variety of magnetic properties [l[ . There are variations from 
itinerant magnetic systems to systems showing characteristics of localized magnetism. In the border between these 
extremes one have the so-called heavy fermions, which show many peculiar and anomalous properties, one of which 
is the co-cxistence of superconductivity and magnetism 0]. One aspect that makes the magnetism of the actinides 
^ ■ unique is the presence of strong spin-orbit coupling (SOC) together with strong exchange interactions for the 5/ 
electrons, which arc the ones responsible for the magnetism. 

From a theoretical point of view, a standard density functional approach, cither in the local density approximation 
(LDA) or generalized gradient approximation (GGA), describes quite well the equilibrium properties of at least 
'■^ ■ metallic systems. However, these functionals are known to underestimate the orbital moments which are induced 
Ph [ by the relatively strong SOC [1, 0, This can be remedied by allowing for the so-called orbital polarization 0], 
responsible for Hund's second rule in atomic physics, either through adding an appropriate orbital depending term to 
the Hamiltonian or by adopting the so-called hDA+U approach In the latter method a screened Hartree-Fock 

(HF) interaction is included among the 5/ states only. 
[ There is a drastic difference between the itinerant magnetism of a 3d shell and that of the 5/ shell. In the former 
' the orbital degrees of freedom are quenched due to the process of hopping between different atoms, while in the latter 
[ case the stronger SOC un-quenches them again. Magnetic ordering is relatively abundant among actinide systems 
• due to the strong exchange interactions, but generally the spin moments are strongly reduced compared to a fully 
[ spin polarized value, which sometimes is ascribed to crystal field effects and other times to hybridization. 
■ This paper will focus on the role of the local screened exchange interactions and it will aim to convincingly argue that 
' these interactions, together with an appreciable SOC interaction, are responsible for the reduced spin polarizations 
as well as for a large orbital moment. This is analyzed in terms of tensor moment polarizations of the 5/ shell. 
Firstly we describe the method we employ, the LDA-ft/ method in its most general form with a minimum of free 
parameters; one. This is accomplished by using screened Slater parameters together with an interpolating optimal 
double-counting (DC) term Q. This last degree of freedom can be chosen to be, for instance, the lowest Slater integral 
U, which is used as a varying parameter. This general LDA-|-?7 scheme has been implemented within the state-of- 
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Ci ■ the-art full-potential augmented plane- wave plus local orbitals (FP-APW-l-Zo) method 13, U] in the Elk code [l3l 
This linearized augmented plane-wave method is one of the most accurate scheme to treat the complicate behavior 
of 5/ electrons in actinides [13, [l^- Secondly, we present an analysis method for the resulting ground state. This 
analysis is based on an exact multipole decomposition of the density matrix as well as of the HF exchange energy. 
We have tested the generality of the method by considering actinides systems with different degree of localization of 
their 5/ electrons; the itinerant features in U compounds, the more localized behavior in Pu compounds and finally 
the intermediate one in Np based materials. The calculations show reasonable comparison with experiments as well 
to other beyond-LDA calculations. By applying our tensor moment decomposition, we observe several clear trends 
regarding the favoured polarizations in these actinide systems. For instance it is clear that the Hund's rules break 
down in the sense that other channels are more significant than the spin polarization. 
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METHOD 
A General Form of LDA+C/ Method 

In the most general version of LDA+U @, 0, EE] a HF correction to the energy enters with the form 

Ehf = ^YPac [{ab\g\cd) - {ab\g\dc)] ptd , (1) 

abed 

where pab is one element of the density matrix with dimension D = 2{2i + 1), which acts as an occupation matrix 
and a, b, c and d are single electron states. The interaction term is of the form 

{ab\g\cd) = / v4(l)V'l(2)g(ri2)Vc(l)V'd(2)d(l)d(2) , (2) 



with one-electron states a with wavcfunction V'a(l) = Riifi)Yemai^i)Xsai^)i where the relevant quantum numbers, 
771 and s, are the magnetic quantum number and spin component, respectively. The interaction can be expanded in 
a series 

oo 

g{ri2) = Ygk{ri,r2)Pk{cos0i2) , (3) 
where the Legendre function P^. in turn can be expanded by the use of the addition theorem for spherical harmonics 

k 
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Pkicos0,2) - Yk,i^i)Yk,in2) . (4) 

q— — k 



Calculation of Slater Integrals through a Yukawa Potential 

The radial part of the interaction is then contained in the Slater integrals 

= / dnrlRj{rMn,r2)RUr2y2dr2 . (5) 



For a / shell there are four independent parameters. F'-^^ = U, F'^'. F^'^^ and F'^^\ It is very unpractical to stay 
with these four parameters. A common practice within LDA+C/ or HF calculations 0, [H, is to have the screened 
Slater parameters determined by the choice of two linear combinations of parameters, U and J, and by fixing two 
ratios, Ai = F^^^/F^^^ and A2 = F^^^/F^^^. In the present work we will instead follow the ideas of Norman 
and calculate the Slater parameters directly from a screened Coulomb interaction in the form of a Yukawa potential 
gii'i,r2) = e~-^''i2/ri2. Then 

i?fe(ri,r2,A) = -(2fc + l)Ajfe(*Ar<)/iW(iAr>) , (6) 

where jk is a spherical Bessel function, /i^^^ is spherical Hankcl function of the first kind and r< and r> are, respectively, 
the smaller and the larger radius entering in the double integral in Eq. ([5]). This type of approach has two advantages; 
it determines the ratio between the different Slater parameters in a more realistic way than by choosing U and J 
individually, and there is only one independent parameter, the screening length A. 

Since in APW+Zo basis set Re{r, e) is energy dependent, we decided to use the energy e at the center of the band 
of the localized shell £. We set the atomic muffin-tin (MT) radius to a value large enough such that the integrals 
in Eq. ([5]) are well converged. In the upper part of Fig. [T] we plot the calculated Slater parameters for US. These 



values are in perfect agreement with ones calculated for the ion by Norman 17[ with the same screened potential. 
In the lower part of Fig. [1] we compare the Slater integrals ratios Ai and A2 obtained for US with the fixed ratios 
commonly used in most LDA-I-C/ studies @, There is a good agreement only at small values of the screening 
length. For A > 0.5 a.u.^^ there start to be a significative difference that turns out to be relevant in the calculation of 
US magnetic moments. If the F^'^'^s are calculated individually with Eq. (O and Eq. ©, the spin moment (Afgp) and 
orbital moment (Morb) show a more dramatic variation as a function of U (or A) than ones determined by fixing Ai 
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and A2 (see Fig. [2])- This fact indicates how relevant the individual determination of every Slater integral might be in 
many systems. Finally we verified that one can reproduces the LSDA moments in the limit of large screening length 
(^(fe) _> for A — > 00), i.e. by slowly changing the parameter A one increases the localization of the 5/ electrons from 
the LSDA limit in a continuos way. 

A parameter-free method to screen the Slater parameter has been suggested by Brooks [l^, where the screening 
parameter A is identified as the Thomas-Fermi screening, which depends on the local charge density. A more appropiate 
but time-consuming scheme to calculate the Slater parameters from the screened Coulomb potential is within the RPA 
approach, as it has been recently accomplished by Solovyev et al. [20l |. 



Calculation of LDA+L^ Potential 

The contribution to the orbital potential from the LDA-f [/ correction is defined as 

Vij = -J— = XI tO'^^lsl*^) - Pab ■ (7) 

ab 

We note that the potential so defined is the complex conjugate of the one sometimes stated in literature. The correct 
definition becomes crucial to evaluate off-diagonal spin terms in all calculations in which the coordinate system is not 
rotated to the local one, i.e. in all calculations in which the density matrix is not diagonal. 



Double-Counting Corrections 

A major obstacle in the LDA+C/ approach is that the electron-electron interaction has already been included in 
LDA potential, thus a simple addition of the orbital dependent HF potential would lead to DC terms. One may want 
to individuate those terms in the LDA potential that correspond to the interaction already considered in the HF 
Hamiltonian and subtract them. A direct connection between the two formalisms is not possible and in addition it 
would not be useful. In fact LDA approximation treats very accurately spatial variations of the exchange-correlation 
potential but it neglects the orbital dependence of the Coulomb interaction. Thus the best recipe would be to identify 
the mean- field part of the HF potential and subtract it, leaving only an orbital dependent correction to the mean- 



field type LDA potential. Czyzyk and Sawatzky 21[ suggested a prescription that is exact in the case of uniform 
occupancies (around-mean-field, AMF) and that would be realistic for weakly correlated systems, however not exact 
because of the presence of the crystalline field. The AMF correction is implemented by redefining a new density 
matrix without the charge n and the magnetization m contributions 

n = Tr/9 (8) 
rfi — Trap , (9) 

in the following way 

Pab = Pab - (Sabn + aab ■ m) / D . (10) 
The AMF-double-counting corrected LDA+C/ energy and potential terms become in our formalism 

E^Tdc = ^ E -^ac [{ab\g\cd) - {ab\g\dc)] pm , (11) 

abed 



yAFM ^ „ ^ [{ja\g\ib) - {ja\g\U)] {S,,n + a,, • m)/D - ^ [{ja\g\zb) ~ {ja\g\bi)] p,, . (12) 

ab ab 

For strongly correlated systems it exists another prescription for the DC, the fuUy-localized-limit (FLL) that 
would correspond to subtract the average effect for a localized state, with integer occupation number. The most 
general expressions for energy and potential are 

E^^^ ^ {2Un{n - 1) - 2 Jn(n/2 - 1) - Jm ■ m} /4 , (13) 
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U{2n - 1) J{n - 1) 



% + ^ . (14) 



Most of LDA+C/ calculations use one of these approaches, while the real occupation numbers lie somewhere between 
the two limits. Petukhov et al. [3| proposed a linear interpolation between these two limits (INT DC), 

= aEl^J^,^ + (1 - a)E§^l , (15) 

V^^^aV^^ + {l-a)V^^ . (16) 

in which the parameter a is a material dependent constant determined in a self-consistent (SC) way following a 
constrained-DFT philosophy. In our formalism the expression of a of Ref. 0] is generalized to take into account the 
off-diagonal spin terms of the density matrix, 

DTvp-p 

" = T, 2 2 ' (1 

Vn — 71^ — 

where n and m arc defined in Eq. ([5]). 

In the present study we prefer to use the INT DC approach for two reasons. Firstly, it reduces one further free 
parameter. The results do depend on the choice of DC, but if wc stay consequently with the INT DC this degree 
of freedom is gone since wc can in principle treat both more itinerant and more localized system. Secondly, it turns 
out that the use of the INT DC is very important to reproduce the correct magnetic structure of monopnictides Pu 
compounds. In Fig. [3]we compare Mgp and Moib calculated for the ferromagnetic PuS [if] with ones calculated for 
the paramagnet PuP [l|. For AMF type of DC, Msp and Aforb of both compounds decrease dramatically until they 
disappear at C/ « 2.0 eV. Instead, by using FLL DC, Mgp and Morb decrease significantly but they never disappear 
for any value of U. Finally, only by using the INT type of DC, we find a range of values for U (J7 > 4.0 eV) in which 
PuS is magnetic and PuP is non-magnetic. 

In conclusion wc have implemented the LDA4-J7 method in the most general form, taking into account the off- 
diagonal spin terms of the density matrix and the correct definition of the potential for those terms. By using the 
interpolated DC of Petukhov et al. Q and by calculating SC the F^'^'^s with a Yukawa potential, our LDA-I-J7 approach 
has only one free parameter left, i.e. the screening length A or if preferable U. 



Multipole Representation of LDA+(7 Energy 



The formalism up to now is standard and have been used several times before, although some studies have neglected 
the spin-mixing terms in Eq. ([T]). Within this formalism the density matrix plays a crucial role. In the following we 
will decompose this fairly large matrix into the most important and physical relevant terms. We will find that this 
decomposition largely simplifies the analysis, as well as it gives many new insights into the magnetism of the actinidcs, 
where the SOC has a crucial contribution. The interaction term in Eq. ([T]) have been studied in detail, for example 
by Slater Racah [13, [3, 0, [2^ Condon and Shortly [27| . By expanding the interaction in spherical harmonics 
and by making use of the Wigner 3j-symbols [27[ the interaction can be expressed as 
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(18) 
(19) 



The spin dependence is given by the two delta functions of the spin quantum numbers of states a, b, c and d. The 
radial dependence is confined in the Slater integrals F^'^^ and the Wigner 3j-symbols take care of the angular part 
of the integral. We now introduce a multipole momentum tensor w*^' defined as the expectation values of a tensor 
operator v'^, 
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{mt,\v';\ma) - (-1)^-'"^ 
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(20) 
(21) 

(22) 
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where the tensor component index x runs from — fc to k. The spin independent part of the HF energy (the Hartree 
term) can be rewritten as 

Eh = ^ 2^n,k I Q Q I -w . (23) 

In order to take care of the spin dependence we may introduce a double tensor 

.Bg = ii-..J(;^ (24) 

where the index y runs from —p to p. It is easy to verify that w'"'" = w*^. Then the exchange energy can also be 
written as a function of the tensor components rather than the density matrix (see Appendix), 

k=0 ki=0 



£ k iV ( £ £ ki^ ^ 



/ \ £ £ k 



^w'^'i'^-w'^if , (26) 

p=0 



where the {...} symbol is the Wigner 6j-symbol. Notice that the Wigner 3j-symbols arc defined such that the 
contribution from odd k vanish, so only Slater parameters of even k are needed. This type of expression was derived 
by Racah [2^ . However, since it was derived for atomic configurations only, it has not been fully realised that it is as 
valid for non-integer occupations. 



The coupling of indices — irreducible spherical tensor 



It is useful to introduce the irreducible spherical tensors w*'''"" from the double tensors w*^^ for two reasons. Firstly, 
the double tensors are not true spherical tensors and, secondly, in the presence of SOC the spin and orbital degrees 
of freedom are not longer decoupled. The three-index tensors w'^^'' is defined through a coupling of the indices of the 
double tensor w'^''. 



kpr 



+P-V 



xay 



k r p 



(27) 



where the index r runs from \k — p\ to |fc +p| and where the normahzation factor n^^^ given, as in Rcf. [28| , by 



abc 



{g-2a)\{g-2b)\{g-2c)\ 



(5 + 1)! 



1/2 



5! 



(5-2a)!!(g-26)!!(g-2c)!! 



(28) 



with g = a + + c. 

These tensor moments have a very nice feature, they are proportional to the moment expansions of the charge (k 
even and p = 0), spin magnetization {k even and p = 1), current {k odd and p = Q) and spin current (fc odd and 
p = 1) densities. For instance, w'"'° gives the total charge, w"^^ gives the spin moment and w^*'^ is proportional to 
the orbital moment. 

The exchange energy Ex of the shell £ in terms of the irreducible spherical tensor moments is now 



Ex - -E^^^^ E E &l)!l^^:i±lK^(-i)^Ha 



2 2 



£ k £ 
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^kipr _ ^kipr 



It is convenient to rewrite Eq. 



in a simplified form, 
i?x = E = E Kk^prw'^'"'- ■ w'^^ 

kipr kipr 



pr 



(29) 



(30) 
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where 

. _ f ^<«Pi±i)!M(^±i),_,,..K.^P„?,. ,31, 

2fc=0 

/ I) 1, n ( fi n 1, > 

(32) 

In Eq. ([5D|) the exchange energy of the sheU i is expressed as a sum of independent terms involving different spherical 
tensors. We will refer to these terms as different exchange channels. 

One of the most important types of polarization is the so-called spin polarization (SP), often referred to as Stoner 
exchange or Hund's first rule, which corresponds to a polarization of channel Oil. Since w*'^^ • w"^^ ~ ™spinj S®^ 
that the SP energy E'sp is given by, 





k 
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EsP = E°l' = A'oiims'pin = -\ + J ) T^spin ■ (33) 



Hence, the so-called Stoner parameter /, defined by Esp = —jlm^p-^, is given by 

/^l^ + J. (34) 

as it is already known. In this expressions we have adopted the convention to use certain linear combinations of the 
Slater parameters, the "Hubbard" -parameters, U = and J, that for d and / electrons is given by, respectively, 

jd ^ j_^pi2) ^ pH)^ ^ (35) 

jf ^ 1f^^) + + J!Li5^(6) . (36) 

45 33 1287 ^ ' 

In our multipolc expansion in Eq. (|30p it is also included an exact formulation of the orbital polarization (OP) 
exchange energy Eqp, 

7T1 _ TTllOl I pllO _ 7,- „,101 ,,,101 I 7,' ,,,110 ,,,110 fn^\ 

il/OP — -c-jif + -l^x — J^ioiw -w +A110W -w , [61) 

where i^m i = 3i^ri 10 and w^^"'^ -w^"^ = "T?^- This expression is a sum of two terms, one that breaks the time-reversal 
symmetry [3l[, 101 (OP-odd), and a second one that does not, 110 (OP-even). Consequently 101 is associated with 
the presence of an orbital moment, while 110 is compatible with a non-magnetic solution [39|. Finally, the prefactor 
Kioi has a simple expression in terms of Racah parameters [isj : this is for d and for / electrons, respectively, 

^101 ^ > (38) 

RESULTS 

In this section we shall apply our method to calculate the magnetic structure of some metallic U, Np and Pu 
compounds, for which the behavior of 5/ electrons varies from itinerant to more localized. 

Application to US 

We selected US as a prototype system to compare the results obtained by calculating the F^^^'s through a screened 
Yukawa potential with other common procedures 0, , for which U and J are provided as inputs and the ratios 
Ai = F^^^/F^"^^ and A2 = F^^l/F^'^^ are fixed. US shows a ferromagnetic order with Tc = 178 K and a strong 
anisotropy along [111] direction 11]. Many experiments have proven 5/ electrons of US to be itinerant 34, Hi,!!!]. 
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Neutron scattering experiments have measured the 5/ total moment of US to be 1.7 /is [37|, while measurements of 
the magnetization in the bulk have shown that the total moment per formula unit is 1.55 hb (38| . On the theoretical 
side there have been many investigations of magnetic properties of US, making this system a good benchmark case 
for our method. Recent LSDA calculation with SOC correction seem to underestimate the orbital moment while 
orbital polarization Q and Hartree-Fock tight-binding calculations reproduce the correct size of Morb We have 
investigated the magnetic and electronic structure of this compound using our optimized LDA+SOC+iJ method and 
its multipole decomposition. As anticipated in the method section, we compared the behavior of magnetic moments 
obtained by calculating F^^'^s through a Yukawa potential (MY"'^^^'^(A) and Mj^l^''^'' {X)) with ones determined by 
fixing the ratios Ai and A2. In the fixed ratios calculations, firstly, we fixed the value of J to 0.46 eV [l3| for any 
values of A (Fig. [2]), secondly, we varied J as function of A to the value J(A) determined by the screened calculation 
(Fig. U]). If the F'^-'^s are determined individually M^^^^^^{X) and Mj^^^^^{X) change dramatically as function of A 
(or U), while the moments are much more constant by fixing Ai, A2 and J. If then J is varied to J(A), and only Ai, 
A2 are kept constant, the moments vary in a similar way to Mj^"'^™(A) and M^^''™^(A). This because in the latter 
case the values of F'^'^^s become fairly similar to ones calculated with the screened potential, as we show in Tab.[Tl 
Let us now analyze the decomposition in multipoles of the HF-Fx energy in Eq. dSO]), that we report in Fig. El We 
found three main contributions: the SP term Oil, the OP-even term 110 and finally the high multipole of m(r) 615. 
For small U the SP Oil is the dominant contribution until it starts to rapidly decrease. At t/ « 0.6 eV the OP-even 
110 crosses Oil and it becomes the largest term. Such a term has been detected in all our LDA-I-SOC-I-C/ calculations 
on actinides and it has been found responsible for the vanishing of magnetic moments in 5-Pu both in a LDA+U 
investigation [sl] as well as in DMFT study [i^. Finally at C/ « 0.6 eV the 615 contribution becomes the largest. 

Few other terms are present but they are much smaller than the others; 101 , OP-odd term, that refers to the Fx 
associated to the presence of an orbital moment, 211 and 413, related to the multipoles of the m(r), 505, corresponding 
to an high multipole of the current density. 

Application to Np compounds 

The 5/ electrons in Np compounds are expected to have intermediate features between the itinerant behavior in U 
compounds and the more localized behavior in Pu compounds. We calculated the electronic and magnetic structure 
of two Np compounds, NpN and NpSb. NpN is reported to be a ferromagnct with easy axis along [111] direction, 
Tc = 87 K and a total magnetic moment of 1.4 /xb [H, NpSb is reported to be an anti-ferromagnet in AFM-3k 
structure and easy axis along [111] direction with a total moment of 2.5 /ib [ll- In Fig. [5] and Fig.[7]we plot Mgp and 
Morb as function of the parameter U and the corresponding screening length A. Both compounds are magnetic for all 
values of A. In the same figures we also report the different multipole components of Fx as defined in Eq. ((30)) . For 
< [/ < 1.5 eV the SP Oil term dominates in both compounds. For U > 1.5 eV, 110 and 615 become the largest 
terms. In both materials there is also a significant OP-odd 101, that is related to the presence of an orbital moment. 
Finally, in NpN wc observe a significant 211 contribution, that corresponds to an high multipole of m(r). We note 
that the relevant channels in these Np compounds are similar to ones present in US, however the 110 term is larger 
and the 615 term is slightly smaller. 

Application to Pu compounds 

We have investigated the magnetic structure of some Pu compounds belonging to the actinides monopnictides, 
whose 5/-electrons are expected to have more localized features compared to Np and U compounds. These materials 
have the fee NaCl crystal structure and the determination of their magnetic structure has been object of many 
experimental studies [ij. PuS, PuSe and PuTe, are paramagnets like S-Pu; PuSb, PuP, PuAs are ferromagnets; 
PuBi is an antiferromagnet. We have found that the use of the INT DC is essential to correctly reproduce the 
magnetic structure of those compounds. In addition we reproduced the non-magnetic ground state of the high-Tc 
superconductor PuCoGas and wc analyzed it through the multipole decomposition of Fx in Eq. ([30]) . Wc already 
discussed in the method section the ability of INT DC to correctly describe the magnetic structure of two prototype 
Pu compounds, the ferromagnct PuS and the paramagnet PuP (see Fig. [3|). The results of a similar calculation for 
the ferromagnets PuAs, PuSb, for the antiferromagnet PuBi and for the paramagnets PuSe and PuTe are summarized 
in Tab. |TT1 By using the INT type of DC, Mgp and Aforb decrease significantly faster for non-magnetic compounds 
until they disappear, while for magnetic compounds the moments decrease slower and do not disappear. 
Let is now discuss the decomposition in multipoles of the HF-Fx- Again we will refer to the two prototype materials 
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PuS and PuP calculated with INT DC. In both compounds the dominant term is SP Oil for < (7 < 1.5 cV. With 
increasing U the OP-cven term 110 starts to increases and it soon becomes the dominant contribution, for U > 1.5 cV. 
In PuS the OP-even 110 term completely takes over the SP Oil contribution, as we already found in the case of 5- 
Pu [3§|. Also in PuP Oil decreases, however it does not become zero and few other significant channels are opened. 
The most relevant one (after 110) is the 615 contribution that is related to an high multipole of m(r). However, 
the 615 term is lower than one present in magnetic U and Np compounds, while 110 channel is clearly larger. This 
mechanism is clearly shown in Fig. [3] and Fig.[8l where we compare the magnetic moment and HF-i^x contributions 
of PuS, non-magnetic, and PuP, ferromagnetic, as function of the parameter U and corresponding screening length A. 
Finally we applied our one-parameter LDA-I-J7 method to the superconducor PuCoGas. Previous LDA+U calculations 
of this compound have stabilized a non-magnetic solution with AMF type of DC, C/=3 cV, J=0.6 eV and fixed Ai 
and A2 4l|. We also stabilized a non-magnetic solution for U~3.2 eV, corresponding to A = 1.79 a.u.~^, (see Fig. [S] 
and Table lU ) . In the Ex multipole decomposition the OP-cvcn 110 term, that corresponds to an enhancement of 
SOC interaction, is again the dominant one taking over the whole SP Oil contribution once the system becomes non- 
magnetic. This mechanism is completel y a nalogous to the one present in the others paramagnetic Pu monopnictides 
we calculated in this work and in 6-Fu 133 . 



CONCLUSIONS 



The purposes of this paper are two-fold. Firstly, we advocates an approach to LDA-I-J7 calculations that reduces 
the number of free parameters in a well-defined way. This involves the use of screened Yukawa potential to describe 
the HF interaction within an atomic shell with a single parameter. In addition we argue to systematically use the 
INT approach to DC, which takes away the ever existing choice between AMF or FLL, especially since the results 
depend on the choice of DC. This approach ought to facilitate fast and systematic LDA+U calculations. The results 
of which would be comparable between different computations, without having to dwell on the used values of e.g. J. 

Secondly, we present a method, the decomposition in tensor moments of the density matrix, that facilitates the 
analyzes of the results from an LDA+C/ calculation. This is an exact approach which gives both the different 
polarization channels as well as the corresponding Hartree and exchange energies. 

Finally, we apply these combined approaches to a few metallic actinide systems whose 5/ electrons show different 
degrees of localization; from the itinerant features in US to the more localized behavior in Pu compounds and the 
intermediate one of Np based materials. The results show that our one-parameter approach catches the ground state 
properties of these system, like spin and orbital moments. Here it is worth mentioning that an even better agreement 
with experiments can be obtained by tweaking the individual Slater parameters differently for the different systems. 
But that would rather reduce than increase the physical understanding of these systems. 

Most importantly, our results signal that these systems are not as well understood as believed. For instance, Hund's 
rules are far away from fulfilled, since the SP docs not play the dominant role as always assumed. On one hand the 
LDA calculations show large SP which leads to an overestimation of the spin moments. On the other hand, when 
including the HF interaction of the LDA-t-C/ functional, another polarization channel dominates, the one of spin- 
orbital currents, the OP-even 110 channel. The tensor moment associated with such channel, w^^°, is even close to 
the saturation value of — Iw*^*"^ for all systems considered. Large values of this tensor moment are in good accordance 
with the rising number of experiments that point to large w^^" for actinide systems 
to a large enhancement of SOC in all these systems by the HF exchange interaction. The SP Oil-channel, however, 
decreases drastically, and in some Pu compounds even disappears. Instead, to some surprise, a third polarization 
plays a large role, the 615-channel; this corresponds to an high multipole of the magnetization density. 



44l . |45| . |46| . This fact corresponds 
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APPENDIX 



The transformation Eq. ([24]) can be reversed through orthogonahty relations of the Wigner 3j-symbols. Hence, the 
density matrix p can be expanded in the double tensors w'^*P', 



^ (40) 

-Sc Ui Sa J -^'^^ 



When this form, together with Eq. pS)) . arc inserted in the exchange part of Eq. ([T]), we get an expression where the 
complications essentially arise from the orbital summations in the factor 

^ V ^"^c x\ ma I \ -ma -q md I \ -mj. X2 mt J \ ~mb q mc J 

maTribnicrridq 

where the indeces on k and x stem from the two different density matrix expansions. The spin dependence of the 
exchange energy is simpler due to the Kroncckcr delta-symbols in Eq. (jlSp . i.e. the relevant factor becomes 



SaSbScSd 



^ V ~Sa ~yi Sb J \ -Sa V2 Sb I ^ ' ' ^^^ ' 



V1-V2 J 



(43) 



by use of orthogonality relations. To simplify Q we start with an identity for the Wigner 6j-symbols, see e.g. Ref. [18, 
and we reshuffle a little, 



e e kiW e e fci 

£ i k ] \ mb —mc —Xi 



nianidq 

£ k £ \f i k i \f £ ki £ ' 
md -q mb J \ -mc q ma J \ ma -Xi -ma 



(45) 



f £ £ A:2 \ 

Then we multiply by ) and sum over mb and mc, by obtaining 

\^ mb —mc X2 J 

£ ^ fci 1 / ^ £ k^ \f £ £ k2\ .^g. 

V ii)ki+t^+'^i+<i ( ^ k, £ \/ £ k e \( e k2 £ \( £ k £ ' 
^ V ""^d ma ) \ -ma -q mc ) \ -mc X2 mb ) \ -mb q m. 

niambmcmdq 

On the LHS we use an orthogonality relation for the Wigner 3j-symbols and on the RHS we identify Q from Eq. ([4T|) 

I \ \ |(2fci + = g. (48) 

The phase factors of Eqs. (|43p and (|48p are then used to form the scalar product of the double tensors in the final 
form of exchange energy, Eq. pSj) as 

^k,p, . ^k,p, ^ Y.^-ir'+y^wl\i\wtT,-y, ■ (49) 
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X (a.u ) 

FIG. 1: Upper panel: Slater parameters i^*^"' (full black line), (red line with squares), 

^(4) ( 

green line with diamonds), 

F^®-* (blue line with circles) and parameter J times 10 (dashed black line) as function of screening length A of the Yukawa 
potential in US. Lower panel: comparison of Slater parameters ratios Ai = f'^'^/F'^' (red lines) and A2 — F'^'/F*-^-* (black 
lines) calculated by a screened Yukawa potential (full line with squares for Ai, full line with circles for A2) with the fixed ones 
of Ref. 0, llll (full hue for Ai, dashed hne for A2) and of Ref. [11, [l3l (dashed-dotted line for Ai, dotted hne for A2). The 
radial wave-function used for the calculation of Slater integrals in Eq. (O corresponds to the energy at the center of 5/ band 
and the MT radius of U is EEt = 2.79 a.u.. 



X (a.u. ^) 




I 1 I 1 I 1 I 1 I 1 I I 

0.4 0.8 1.2 1.6 2 

U (eV) 

FIG. 2: Spin (upper panel) and orbital magnetic moment (lower panel) of US, a = 10.36 a.u., calculated with LDA+SOC+(/ 
approach and INT DC (see DC corrections in the method section). We compare the results obtained by calculating the F'*'s 
with a screened Yukawa potential (full black line, screening length in the upper axis) with ones obtained by fixing the ratios 
Ai and A2 to ones of Ref. [3, [3l (dashed red line) and to ones of Ref. [l^ . \it\ (dashed-dotted blue line) . In the fixed ratios 
calculations we fixed the parameter J=0.46 eV to one of Ref. [T^ . 

The MT radii of U and S are set, respectively, to Rmt = 2.79 a.u. and i?MT = 2.13 a.u.. The parameter R^rj,\G + fc|max, 
governing the number of plane waves in the APW+lo method, is chosen to be 9.76. The BZ is sampled with 1728 fc-points 
uniformly spaced. 
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12 3 4 

U(eV) 

FIG. 3; Spin Msp (red lines) and orbital A/orb (blue lines) magnetic moments of paramagnetic PuS (full line with squares 
for Msp, full line with circles for Morb), a = 10.46 a.u., and ferromagnetic PuP (dashed line for M^p, dashed-dotted line for 
Morb), a — 10.49 a.u., as function of U for different types of DC. The Slater parameters are calculated SC by using a Yukawa 
potential with screening length A reported in the upper axis. For AMF DC the moments vanish in both compounds at (7 ~ 2.0 
eV, for FLL DC both compounds stay magnetic for all values of U. Only for the INT DC there is a range of U (C/ > 4.0 eV) 
for which we obtain the experimental magnetic structure of both compounds PuS becomes non-magnetic while PuP stays 
magnetic. The green lines refer to the 5/ charge of PuS (full line with triangles) and PuP (dotted line). The black lines refer 
to INT DC factor a of PuS (full line) and PuP (dashed-two-dots line). The MT radii of Pu and S(P) are set, respectively, to 
~ 2.82 a.u. and = 2.15 a.u.; the basis set cut-off -R^tIC + ^]max is set to be 9.1. The integration in the BZ is 

performed with 1728 with fc-points uniformly spaced. 
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U (eV) 

FIG. 4: Spin (upper panel) and orbital magnetic moment (lower panel) of US, a = 10.36 a.u., calculated with LDA+SOC+(/ 
approach. We compare the results obtained by calculating the F^'^^'s with a screened Yukawa potential (full black line, screening 
lenght in the upper axis) with ones obtained by fixing the ratios Ai = F'^' /F^^'> and A2 = i^'®' /F^^' to ones of Ref. i, M 
(dashed red line) and to ones of Ref. [l^. [itI (dashed-dotted blue line). In the fixed ratios calculations we varied J to the one 
determined by Yukawa potential at every different screening length, J(A). Other calculation details are the same as the ones 
reported in Fig. [21 



X (a.u. ^) 




U(eV) 

FIG. 5: Relevant exchange energy channels in Eq. ((SO]), -Ex of US calculated with LDA+SOC+(7 method and INT DC. 
The Slater parameters are calculated SC by using a screened Yukawa potential with screening length A reported in the upper 
axis. Other calculation details are reported in Fig. [2] 
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FIG. 6: Upper panel: spin (red line with squares) and orbital (blue line with circles) moment of ferromagnetic NpN, a = 9.25 a.u.. 



Lower panel: relevant exchange energy channels in Eq. E^^'"' , of NpN calculated with LDA-|-S0C-|-(7 approach and INT 

DC. The Slater parameters are calculated SC by using a screened Yukawa potential with screening length A reported in the 
e MT radii of Np and N are set, respectively, to = 2.60 a.u. and -Rmt — 1-80 a.u.; the parameter 

is set to 9.4. The integration in the BZ is performed with 1728 fc-points uniformly spaced. 



upper axis. 
-Rmt \G -\- k 




U(eV) 



FIG. 7: Upper panel: spin (red line with squares) and orbital (blue line with circles) moment of NpSb in AFM-3k structure, 
a = 11.82 a.u.. Lower panel: relevant exchange energy channels in Eq. (|30|l . E^^^ , of NpSb calculated with LDA-I-SOC-I-C/ 
and INT DC. The Slater parameters are calculated SC by using a screened Yukawa potential with screening length A reported 
in the upper axis. The MT radii of Np and Sb are set, respectively, to B!^r^ = 2.95 a.u. and -Rmx ~ 2.65 a.u.; the parameter 
^mtI^ + felmax is Set to bc 9.2. The integration in the BZ is performed with 1728 fc-points uniformly spaced. 
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FIG. 8: Relevant exchange energy channels in Eq. (|30p . E^^^ , for paramagnetic PuS, a — 10.46 a.u., and ferromagnetic PuP, 
a = 10.49 a.u., as function of U for INT DC. The Slater parameters are calculated SC by using a Yukawa potential with 
screening lenght A reported in the upper axis. Other calculation details are reported in Fig. |3] 
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FIG. 9; Upper panel: spin (dashed red line) and orbital (dashed-dotted blue line) magnetic moments of high- Tc superconductor 
PuCoGas, a = 7.84 a.u and c = 12.57 a.u., calculated with LDA+SOC+L'' and INT DC, as function of parameter U. The 
interpolation factor a is plotted with a dotted black line. The Slater parameters are calculated SC by using a Yukawa potential 
with screening lenght A reported in the upper axis. We used an AFM-lk structure with q = (0, 0, 1/2). Lower panel; relevant 
exchange energy channels in Eq. (|30p . E-^'"" , of PuCoGas. The muffin-tin (MT) radii of Pu, Co and Ga are set, respectively, 

=2.2 a.u.. The parameter R^^\G + fc|max is set to 9.5. The integration in the BZ is performed 



to m 



2.7 a.u. and R 



,Co,Ga 
■MT 



with 576 fc-points uniformly spaced. 
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TABLE 11: Table of screening length A used to determine Slater integrals through a screened Yukawa potential in LDA+S0C+(7 
calculations of Pu compounds with INT type of DC. We report the smallest value of U and corresponding screening length 
A that are necessary to obtain vanishing magnetic moments (less than 0.1 fj,B for both Msp and Morb) in non-magnetic Pu 
compounds. For the first three non-magnetic compounds PuX, with X=S,Se,Te, we considered the corresponding ferromagnetic 
compounds PuY, with Y=P,As,Sb, where X and Y are chosen to be neighbours elements in the periodic table. We then report 
the magnetic moments of PuY using the same value of U that makes the corresponding compound PuX non-magnetic. For the 
antiferromagnet PuBi we used the same value of U of PuSb. We also write down the INT DC interpolation factor a that is 
determined SC. Values of experimental magnetic moments are from Ref. [Jj. 
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[a.u-^] 
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PuS 
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1.74 


0.75 










PuSe 


3.6 




1.57 


0.79 










PuTe 


4.1 




1.38 


0.83 










5-Pn 
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0.78 










PuCoGas 
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1.79 


0.69 
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] Morb [ms] 


Mtot [fJ-B 


1 M-P [/is] 


Q 


PuP 
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1.60 


-1.52 


0.08 


0.75 


0.82 


PuAs 
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-1.54 


0.13 


0.64 


0.86 


PuSb 
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1.40 


1.57 


-1.47 


0.1 


0.67 


0.9 
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1.38 


-1.29 
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